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ABSTRACT 

We have investigated the relationship between the X-ray variability amplitude 
and black hole mass for a sample of 46 radio-quiet active galactic nuclei observed 
by ASCA. Thirty-three of the objects in our sample exhibited significant variability 
over a time-scale of ^40 ks. We determined the normalised excess variance in the 
2-10 keV light curves of these objects and found a significant anti-correlation between 
excess variance and black hole mass. Unlike most previous studies, we have quantified 
the variability using nearly the same time-scale for all objects. Moreover, we provide 
a prescription for estimating the uncertainties in variance which accounts both for 
measurement uncertainties and for the stochastic nature of the variability. We also 
present an analytical method to predict the excess variance from a model power spec- 
trum accounting for binning, sampling and windowing effects. Using this, we modelled 
the variance-mass relation assuming all objects have a universal twice-broken power 
spectrum, with the position of the breaks being dependent on mass. This accounts for 
the general form of the variance-mass relationship but is formally a poor fit and there 
is considerable scatter. We investigated this scatter as a function of the X-ray photon 
index, luminosity and Eddington ratio. After accounting for the primary dependence 
of excess variance on mass, we find no significant correlation with either luminosity 
or X-ray spectral slope. We do find an anti- correlation between excess variance and 
the Eddington ratio, although this relation might be an artifact owing to the uncer- 
tainties in the mass measurements. It remains to be established that enhanced X-ray 
variability is a property of objects with steep X-ray slopes or large Eddington ratios. 
Narrow-line Seyfert 1 galaxies, in particular, are consistent with being more variable 
than their broad line counterparts solely because they tend to have smaller masses. 
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1 INTRODUCTION 

Variability was discovered in the X-ray emission from ac- 
tive galactic nuclei (AG Ns) roughly three decades ago 
fe.g.. lMarshall et alll98ll and references therein). EXOSAT 
subsequently obtained well-sampled light curves on time- 
scales of minutes to days, and the power spectra gener- 
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ated from these light curves were described as a power- 
law P oc v~ a with a steep 'red-noise' index of a ~ 1.5 
and an amplitude inversely proportional to the l uminos- 
ity llLawrence fc Papadakislll993l : iGreen et al]ll993lh It was 
clear that this power-law must break at some lower fre- 
quency, or the power would diverge, and some evidence 
for this was found using longer-term archival observations 
jMcHardvlll988l : iPapadakis fe McHardvlll99rJ) . It was not 
until the launch of R XTE, however, that thi s break was 
measured definitively jEdelson fc Nandralll999l) . 

A number of high quality power spectra have now 
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been obtained, pri marily using RXTE and XMM- 
Newt on data (e.g.. lUttlev et alJ l2002t IMarkowitz et all 
\ 200A IVaughan et all l2003bl: iMcHardv et all 120041: 
lUttlev fc McHardvl l20plh . These have shown breaks to 
be common and emphasized the similarity of AGN power 
spectra to that of the black hole binary Cyg X-l. In the 
low/hard state, the power spectrum of Cyg X-l exhibits a 
twice-broken power-law which breaks from a slope of a ~ 
to 1 at the 'low- frequency break' (vlfb) and from a ~ 1 to 
2 at the 'high-frequency b reak' (z^hfb), with ^hfb ~l-6 Hz 
(e.g.. lBelloni fc Hasi ngcr 199oQ). In the high/soft state, the 
power spectrum exhibits only a high-frequency break, with 
i/HFB ~10-15 Hz (e.g., ICui et alNl997l : iRevnivtsev et all 
2000). Though still a subject of debate, the emerging 
consensus is that we usually see the high-frequency break 
in the AGN power-spectra, although two breaks are appar- 
ently seen in two objects (viz, AKN 564 an d NGC 3783; 
IPapadakis et al.ll2002t IMarkowitz et al.ll2003ft . In another, 
the narrow-line Seyfert 1 (NLS1) NGC 4051, there is 
no low-frequency tur nover to a = dow n to very low 
frequencies, which led lMcHardv et all (|2004) to hypothesise 
that this object (and possibly all NLSls) resembled Cyg 
X-l in the high/soft state. 

Determining accurate power-spectra for AGN is diffi- 
cult, as it requires high quality data with near-even sam- 
pling. Such data are available only for a limited number of 
objects and are very costly to obtain in terms of observing 
time. It is nonetheless very useful to quantify the X-ray vari- 
ability of AGN to compare with other properties, and nor- 
malised exce ss variance, denoted as ct^xsj i s much simpler 
to calculate iNandra et alJll997al) . An anti-correlation was 
found between excess variance and luminosity for a sam- 
ple of AGNs observed by the Advanced Satellite for Cos- 
mology and Astrophysics (ASCA), confirming the EXOSAT 
results but with a larger sample of objects (iNandra et alJ 
Il997ah . Later work also using ASCA data revealed that, for 
a given luminosity, the X-ray light curves of NLSls exhibit 
a larger excess variance than the c lassical Seyfert 1 galaxies 
iTurner et alJll999t lLeighlvlll999ah . 

iLu fc Yul <200ll) and iBian fc Zhaol i2003l) . again using 
ASCA data, studied the relationship between the excess 
variance (on a time-scale of roughly 1 d) and the black hole 
mass. Those studies revealed an anti-correlation between 
°nxs ari d mass, which is suggestive that this is the pri- 
mary relationship rather than with luminosity. The NLSls 
appeared to follow the same relationship as the other AGN. 

IPapad akis ( 2004:) investigated the relationship between 
excess variance and black hole mass on much longer time- 
scales ( ~ 300 d) using RXTE data on a sample of 10 AGNs. 
The classical Seyfert 1 galaxies followed a variance-mass 
relation that is consistent with a universal power-spectral 
shape as described above for the l ow/hard state of Cyg X- 
1. In the universal model used bv iPapadakisI <j2004l. i^hfb 
is inversely proportional to black hole mass, and the am- 
plitude, when represented in power x frequency space, is 
assumed to be constant. In agreem ent with th e powe r spec- 
trum analysis of McHardy (2004), Papadakis (2004) found 
that the NLS1 NGC 4051 did not follow the same variance- 
mass relationship described by the classical Seyfert Is. The 
excess variance of NGC 4051 was consistent with a singly- 
broken power-law, breaking from a = 1 to 2, with a break- 



frequency 20 times higher than that deduced for the other 
the Seyfert Is. 

These works show that excess variance can be a use- 
ful complement to full-blown power-spectral analysis, and 
have the advantage that they can be applied to a larger 
nu mber, and wider variet y of objects. As has been shown 
bv lVaughan et all (l2003all some caution must be exercised 
when interpreting excess variance measurements, primarily 
due to the red-noise shape of the power spectra and the 
stochastic nature of the variability. Such effects have not 
been accounted for in the majority of previous works. The 
intention of the work presented here is to investigate the re- 
lationship between excess variance and mass in a large sam- 
ple of AGN, improving on these previous studies by fully ac- 
counting for measurement uncertainties, sampling and red- 
noise effects in the calculation of the excess variance and its 
uncertainty. 



2 THE TARTARUS DATABASE AND THE 
AGN SAMPLE 

The Tartarus 1 database contains products f or ASCA ob- 
serva tions with targets designated as AGN (iTurner et alJ 
2001). We selected radio-quiet objects that have data in the 
Tartarus (Version 3) database and also for which we could 
conveniently obtain a measurement of the black hole mass, 
M m . Seyfert 2 objects were excluded from our sample, with 
the exception of NGC 5506 because for this object we are 
confident of seeing the X-ray emission directly jBlanco et ID 
1990). This initial sample comprised 68 AGN. We utilised 
the Tartarus analysis pipeline to extract light curves for the 
objects in this sample. As we describe in detail in the follow- 
ing Section, not all light curves were suitable for our anal- 
ysis. Having screened the available data, there remained 46 
objects for which we could suitably characterise the X-ray 
variability. These objects are listed in Table Note that, 
while a flux limit was not formally applied to our sample, 
the effect of the screening process was to exclude objects 
having a low counting rate. 

Recent progess in measuring black hole masses has 
made possible the work we present here. We preferen- 
tially used the rever beration-mapping mass estimate from 
iPeterson etafl J2004I) . If this was not available then we used 
the mass estimat e as determined from either the stellar ve- 
locity dispersion dGebhardt et alJ jgOOO'l or the empirical re- 
lationship b etween the broad-lin e region radius and 5100 A 
luminosity (IWandel et al.lll999l) . The masses are given in 
Table Q where we also list the method used to determine 
the mass and the corresponding reference. The masses for 
most objects were available in the literature. For 8 objects 
in Table U we obta ined optical spectral infor mation from 
iGjMa pe et alJ i2004l) and utilise d eqn. 6 from iK aspi et all 
fcOOd) and eons. 1 and 2 from IWoo fc Urrvl feOOSft to de- 
termine M m . 

The 2-10 keV luminosity L2-10 kev and hard- X-ray (ei- 
ther 2-10 keV or 3-10 keV) photon index T are also listed 
in Table for those objects in which we detected variabil- 
ity. The majority of L2-10 kev and T values were taken 

1 http:/ /astro. ic.ac.uk/Rosearch/Tartarus 
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Table 1. X-ray spectral and variability information for objects having at least 1 valid light curve segment. The 2-10 keV luminosity 
and hard- X-ray photon index are given for objects in which variability was detected. 



Name 


M m 




r 


Num. 
Seq. 


Num. 

See 




2 

"nxs 
± Boot. Unc. 




lot? (T 2 T ^~ 

u s NXS 
± Total Unc. 


Aloe rr 2 „,~ 
& NXS 

± Total Unc. 


Refs. 




(2) 


(3) 


(4) 


(5) 


(61 




(7) 




(8) 


(9) 


(10) 


MRK 335 


7.15 


43.07 


1.87 


1 


i 


(3 


12 + 1.87) x 10 


-3 


-2.51 + 0.41 


-0.44 + 0.41 


R,l,2 


PG 0026+129 


8.59 


44.53 


1.96 


1 


3 


(1 


31 ± 1.92) x 10 


-3 


-2.88 + 0.66 


0.61 + 0.66 


R,l 


TON S180 


7.09 


43.58 


2.43 


2 


26 


(1 


59 + 0.10) x 10 


-2 


-1.80 + 0.07 


0.21 + 0.07 


L,3,4 


I Zw 1 


7.20 


43.35 


2.40 


1 


1 


(1 


88 + 0.92) x 10 


-2 


-1.73 + 0.39 


0.39 + 0.39 


L,3,4 


F 9 


8.41 


43.91 


1.91 


8 


6 


(3 


49 ± 5.52) x 10 


-4 


-3.46 + 0.70 


-0.16 + 0.70 


R,l,2 


RX J0152.4-2319 


7.87 






1 


2 




< 6.5 x 10~ 3 




< -1.94 




L,5 


MRK 0586 


7.86 


44.07 


2.22 


1 


3 


(2 


57 + 0.75) x 10 


-2 


-1.59 + 0.22 


1.17 + 0.22 


L,6,4 


MRK 1040 


7.64 


42.40 


1.69 


1 


1 


(1 


20 + 0.65) x 10 


-2 


-1.92 + 0.40 


0.62 + 0.40 


S,6,7 


NGC 985 


8.05 


43.50 


1.73 


1 


2 


(3 


47+ 1.76) x 10 


-3 


-2.46 + 0.32 


0.48 + 0.32 


L,5 


1H 0419-577 


8.58 






2 


3 




< 4.31 x 10~ 3 




< -2.12 




L,3 


F 303 


6.37 


43.03 


1.92 


1 


1 


(6 


72 ± 6.03) x 10 


-3 


-2.17 + 0.44 


-0.74 + 0.44 


L,5 


AKN 120 


8.18 


43.88 


1.93 


1 


2 


(3 


78 ± 7.67) x 10 


-4 


-3.42 + 0.91 


-0.35 + 0.91 


R,l,8 


PG 0804+761 


8.84 






1 


2 




< 3.37 x 10~ 3 




< -2.23 




R,l 


PG 0844+349 


7.97 






1 


2 




< 1.17 x 10~ 2 




< -1.69 




R,l 


MRK 110 


7.40 






1 


1 




< 1.63 x 10~ 3 




< -2.55 




R,l 


PG 0953+415 


8.44 






1 


2 




< 8.18 x 10~ 3 




< -1.85 




R,l 


NGC 3227 


7.63 


41.66 


1.52 


2 


4 


(2 


41 + 0.20) x 10 


-2 


-1.62 + 0.16 


0.91 + 0.16 


R,l,2 


MRK 142 


6.76 


43.17 


2.12 


2 


1 


(1 


54 + 1.33) x 10 


-2 


-1.34 + 0.34 


0.37 + 0.34 


L,5,4 


HE 1029-1401 


9.08 


44.44 


1.83 


1 


2 


(1 


02 + 1.21) x 10 


-3 


-2.99 + 0.56 


0.98 + 0.56 


L,6,9 


NGC 3516 


7.63 


43.08 


1.83 


5 


18 


(3 


70 ± 0.45) x 10 


-3 


-2.43 + 0.10 


0.10 + 0.10 


R,l,2 


PG 1116+215 


8.21 






1 


1 




< 1.06 x 10" 2 




< -1.73 




L,6 


EXO 1128.1+6908 


7.02 






1 


1 




< 1.78 x 10~ 2 




< -1.51 




L,5 


NGC 3783 


7.47 


42.90 


1.70 


9 


8 


( 3 


91 + 0.51) x 10 


-3 


-2.41 + 0.13 


-0.03 + 0.13 


R,l,2 


NGC 4051 


6.28 


41.21 


1.92 


2 


6 


(8 


62 + 0.66) x 10 


-2 


-1.06 + 0.09 


0.31 + 0.09 


R,l,2 


NGC 4151 


7.12 


42.62 


1.53 


13 


29 


(2 


79 ± 0.22) x 10 


-3 


-2.55 + 0.07 


-0.51 + 0.07 


R,l,2 


PG 1211+143 


8.16 






1 


1 




< 2.39 x 10~ 2 




< -1.38 




R,l 


MRK 766 


6.54 


42.73 


2.16 


1 


2 


(4 


02 ± 0.48) x 10 


-2 


-1.40 + 0.16 


0.15 + 0.16 


L,6,2 


NGC 4395 


4.11 


39.99 


1.7 


5 


6 


(1 


13 + 0.14) x 10 


-1 


-0.95 + 0.10 


0.17 + 0.10 


L, 10,11 


NGC 4593 


6.73 


42.98 


1.81 


2 


1 


(1 


42 + 0.21) x 10 


-2 


-1.85 + 0.33 


-0.16 + 0.33 


R,l,8 


WAS 61 


6.66 






1 


1 




< 6.95 x 10~ 3 




< -1.92 




L,5 


PG 1244+026 


6.07 


43.03 


2.46 


1 


2 


(2 


60 + 0.62) x 10 


-2 


-1.59 + 0.18 


-0.31 + 0.18 


L,5,12 


MCG-6-30-15 


6.19 


42.72 


2.00 


6 


48 


(4 


16 + 0.13) x 10 


-2 


-1.38 + 0.03 


-0.05 + 0.03 


L,3,2 


IC 4329A 


7.00 


43.59 


1.71 


5 


6 


(2 


36 ± 2.44) x 10 


-4 


-3.63 + 0.47 


-1.70 + 0.47 


R,l,2 


MRK 279 


7.54 


43.66 


1.99 


1 


1 


(2 


32 + 0.84) x 10 


-3 


-2.63 + 0.36 


-0.19 + 0.36 


R,l 


NGC 5506 


7.94 


42.73 


2.08 


1 


2 


(1 


06 + 0.14) x 10 


-2 


-1.97 + 0.23 


0.87 + 0.23 


S,13,8 


NGC 5548 


7.83 


43.41 


1.79 


11 


16 


(9 


42 ± 2.67) x 10 


-4 


-3.03 + 0.14 


-0.30 + 0.14 


R,l,2 


MRK 1383 


9.11 






1 


1 




< 6.33 X 10 ° 




< —1.96 




R,l 


MRK 478 


7.34 


43.50 


2.06 


1 


2 


(6 


14 + 3.75) x 10 


-3 


-2.21 + 0.35 


0.04 + 0.35 


L,3,4 


MRK 841 


8.10 


43.54 


2.00 


3 


5 


(1 


14 + 0.93) x 10 


-3 


-2.94 + 0.38 


0.05 + 0.38 


L,6,2 


MRK 290 


7.05 


43.22 


1.77 


1 


2 


(4 


11 ± 2.15) x 10 


-3 


-2.39 + 0.32 


-0.41 + 0.32 


L,3,7 


IRAS 17020+4544 


6.77 


43.73 


2.37 


1 


2 


(5 


47 + 2.00) x 10 


-3 


-2.26 + 0.28 


-0.54 + 0.28 


L,14,4 


MRK 509 


8.16 


44.03 


1.82 


11 


2 


(5 


75 + 7.17) x 10 


-4 


-3.24 + 0.59 


-0.19 + 0.59 


R,l,2 


AKN 564 


6.06 


43.38 


2.58 


13 


70 


(5 


34 + 0.14) x 10 


-2 


-1.27 + 0.03 


0.00 ± 0.03 


L,3,4 


RX J2248.6-5109 


7.67 






1 


1 




< 1.08 x 10" 2 




< -1.73 




L,5 


NGC 7469 


7.09 


43.25 


1.84 


3 


2 


(4 


68 + 1.60) x 10 


-3 


-2.33 + 0.27 


-0.32 + 0.27 


R,l,2 


MCG-2-58-22 


8.54 






2 


1 




< 1.53 x 10" 3 




< -2.58 




L,3 


objects are listed in 


order of R.A. (1) Object 


name. 


(2) Log 


of black hole mass in 


units of Mq . (3) Log of 2-10 keV luminosity in 



units of erg s — 1 . (4) Hard-X-ray photon index. (5) Number of available ASCA observing sequences. (6) Number of usable light curve 
segments. (7) Mean normalised excess variance with the uncertainty or upper limit as determined from the bootstrap simulations. 
(8) Log of the mean normalised excess variance with the uncertainty as determined by combining the bootstrap uncertainty and 
the derived red-noise scatter. (9) Residuals from the best-fitting universal model with the uncertainty as determined by combining 
the bootstrap uncertainty and the derived red-noise scatter. (10) Method used to determine the black hole mass and references for 
the mass and X-ray spectral properties. The methods, in order of preference, are as follows: R, reverber ation mapping; S, ste llar 
velocity dispersion; L , r elationship between bro ad-line region ra di us and optical lumin osi ty. References : 1. [P eterson et alJ i2004j): 2, 
Nandra et ai]il997bl): 3.lBian fe Zhaol i2003l): 4.lLeighIvl (l999rl): 5.lGrupe et al.U2004|): 6.IWoo fe UrrvH2()02l) : 7.lRevnoldsl jl997tl: 8 . 
Nandra fc Pounds! ill 9941) : 9. iReeves fc Turnerl <2000D : lO. lFilippenko fe HoH2003t) : ll. llwasawa et alJ i200cT ri2. lGeorge et aUfcOOO) : 
13. |Papadakis! <2004l) : ll lWang fc Lul <200j) . The reference for the black hole mass of each object is listed first. 



© 2004 RAS, MNRAS 000.IT1IT31 



4 P. M. O'Neill, K. Nandra, I. E. Papadakis and T. J. Turner 



from iNandra et alJ (Il997rl. iLeiriilvl ill999bll. iRevnoldsl 
Jl997l).lNandra fc Poundsl lll994l). iReeves fc Turned fcOOCl) . 
Ilwasawa et, alJ l)2000h . and lOeorge et alJ (j200Ct) ~ For the ob- 
jects PG 0026+129, NGC 985, F 303, and MRK 279, we 
fitted the available ASCA data to obtain L2-10 kcv an d 
T. The SISO, SIS1, GIS2, and GIS3 spectra were fitted 
simultaneously in the 2-10 keV rest-frame energy range. 
We used an absorbed power-law, with Nh constrained to 
be greater the galactic value which we obtained using the 
NASA HEASARC 'nH' tool. 2 For our luminosity calcula- 
tions we obtained redshifts from the NASA/IPAC Extra- 
galactic Database 3 and used Ho = 75 km s _1 and go = 0.5. 
All L2-10 kcv values collected from the literature were trans- 
formed to this cosmology as required. 



3 EXCESS VARIANCE ANALYSIS 

The number of ASCA observing sequences available for 
each object is shown in Tabled We extracted a 2-10 keV 
combined SIS0+SIS1+GIS2+GIS3 light curve from each se- 
quence. These initial light curves had a resolution of 16 s 
and each bin was required to be fully exposed. The light 
curves were then rebinned to a resolution of 256 s. 



3.1 Excess variance calculation 

For a red-noise process, the variance in a light curve depends 
both on the power spectrum of the variations and also on the 
time resolution and duration of the light curve. This means 
that different On XS measurements are only strictly compa- 
rable if the durations of the light curves are equal. There- 
fore, we sub-divided the light curve from each sequence into 
many segments of similar duration. The advantage of using 
long durations is that the amplitude of variability increases, 
and the number of points used to calculate cr^xs i s a l so 
larger, reducing the measurement uncertainty. On the other 
hand, using short durations has the advantage that more 
light curve segments can be included. We chose a nominal 
segment length of 40 ks for our analysis as a tradeoff be- 
tween these considerations. In reality we chose a duration of 
39936 ks which is an integer multiple of our 256 s time bin. 

The sub-dividing of the light curves proceeded as fol- 
lows. First, the earliest 40 ks segment of the light curve for 
a certain observing sequence was selected. Then, beginning 
with the next exposed bin following the end of this first 
segment, another 40 ks segment was selected. This contin- 
ued until the light curve had been completely sub-divided. 
Note that the actual duration of these light curves, which 
we define as the time between the first and last exposed bin, 
can be less than 40 ks because the dividing point between 
segments can occur when there is a gap in the data train. 
We accepted all resulting light curve segments that had a 
duration >30 ks. 

To ensure Gaussian statistics, we required each 256 s 
bin to contain at least 20 counts. The number of counts in a 
certain 256 s bin depends both on the source counting rate 
and the fractional exposure of that bin. We do not wish to 

2 http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 

3 http:/ /nedwww. ipac.caltech.edu/ 



reject entire observing sequences simply because some of the 
bins have a low exposure, but those in which fully-exposed 
bins have < 20 counts should be excluded. If we were to 
remove bins simply because they had a low counting rate, 
we would be biased against observing objects when their 
intensity is low. Selecting according to fractional exposure, 
on the other hand, can remove bins having too few counts 
without introducing this bias, as fractional exposure is not 
related to the intensity of the source. If, for example, there 
are fully exposed bins with < 20 counts, the entire sequence 
was discarded. This is the case when we are dealing with a 
weak source. For brighter sources where only underexposed 
bins have < 20 counts, we excise all bins below some min- 
imal fractional exposure. This gets rid of the non-Gaussian 
bins, allowing us to keep the remainder of the light curve for 
further analysis, but introduces no bias against those times 
when the source is weak due to true flux variability. 

Finally, we further required the truncated and screened 
light curve segments to have at least 20 bins, so that the 
variance could be determined accurately. 

This procedure resulted in 46 objects having at least 1 
valid light curve segment, and a total of 305 valid segments 
in all. The number of segments for each object is given in 
Table The mean durations of the light curve segments 
for each object were in the range 35-40 ks in the observers 
frame. The 48 objects in our sample have redshifts in the 
range 0.001-0.234. Taking into account the redshift of each 
object, the rest-frame mean durations were in the range 30- 
40 ks. We expect the effect of these slightly different dura- 
tions to be small. For a power spectrum with a power-law 
slope of a = 2, the worst case we expect, a ~25 per cent 
reduction in the light curve duration (i.e., from 40 ks to 
30 ks) results in a reduction in the 0n XS of only ~0.1 dex. 
As presented later in this Section, the uncertainties in most 
of our observed a^xs values are a few to several times larger 
than 0.1 dex. Therefore, the ~25 per cent difference between 
the shortest and longest mean light curve duration can be 
neglected and allows us to use more data than would have 
been available if we had imposed a strict limit on duration. 

We tested to see which objects exhibited significant vari- 
ability by performing a chi-squared test. The % 2 correspond- 
ing to the hypothesis of a constant counting rate was deter- 
mined for each ~40 ks light curve segment. Then, for each 
object, we summed all of the y^s and degrees-of-freedom 
(DOFs) to test whether that object is variable. We detected 
variability in 33 objects at the 95 per cent confidence level. 
We then calculated the excess variance in each light curve 
segment with the following expression: 

v 

ffNXS = Jj—^l(. X i ~ ~ a i] (!) 
i — 1 

where N is the number of bins in the segment, Xi and en 
are the counting rates and uncertainties, respectively, in each 
bin, and p is the unweighted arithmetic mean of the count- 
ing rates. For objects with more than one valid segment, the 
unweighted average excess variance was determined. A ma- 
jor advantage of our work is that, given the large number 
of light curves available, there is often more than one valid 
segment per object (see Table 0. Taking the mean o^xs 
of these multiple segments reduces the potentially large un- 
certainty owing to the stochastic nature of the variability 
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(see below). When calculating the mean excess variance we 
used all valid light curve segments for a particular object, 
including those segments that did not, in themselves, exhibit 
variability based on the x 2 test. 

3.2 Estimating the uncertainties in a^xs 

Estimating the uncertainty for excess variance is somewhat 
complicated. A nalytical presc r iptions have been given in 
the literature bv lNandra et alj ll997a|, their correct for mula 
is given by Tu r ner et al. 19991. lEdelson et alJ (l2002h and 
IVauehan et alJ i2003al) . The latter authors also discussed 
the uncertainties in cr|[ XS on the basis of simulated red-noise 
light curves. These uncertainties depend both on measure- 
ment uncertainties (e.g., Poisson noise) in the light curve 
data, and the stochastic nature of the variability: any given 
light curve segment represents just one realisation of a ran- 
dom process, and thus can exhibit a different mean and vari- 
ance from the true value, or another random segment. This 
'noise' uncertainty can be very large, especially for a single 
realisation. One must, therefore, account for this uncertainty 
before apparent differences i n cr^ xs , either in a given source 
l|Nandra fc Pa padakisll200 l]) or in comparing sources (e.g., 
lTui*n^r"e'tHuT^)99l) can be considered robust. The measure- 
ment and noise uncertainties on cr^xs are unrelated, so can 
and must be estimated separately. 

To estimate the uncertainty in Onxs owing to mea- 
surement uncertainti es, we used bootst rap simulations (the 
reader is directed to IPress et al] 1200 ll for a discussion on 
bootstrap simulations). Suppose that the observed light 
curve contains TV bins. This light curve is a distribution of N 
counting rates and corresponding Poisson-noise uncertain- 
ties from which we calculate cr^xs- Note that calculating 
°nxs does not depend on the bins being in time-order. A 
bootstrap simulation involves randomly selecting, from that 
distribution, a new set of TV bins. The duplication of bins 
is permitted during the selection process. This, then, results 
in the creation of a slightly different distribution of counting 
rates, and cr^xs can De determined for this new distribution. 
If one repeats the entire process many times, the resulting 
distribution of simulated o"nxs values provides an estimate 
of the uncertainty in o^xs- 

We performed a series of 10000 bootstrap simulations 
to determine the uncertainty in the mean observed o"nxs 
for each object in our sample. Each of these simulations in- 
volved: simulating a new 'light curve' from each valid light 
curve segment, determining cr^xs f° r those simulated light 
curves, and then determing the mean of these simulated 
°nxs values. We were thus able to generate 10000 simu- 
lated values of the mean cr^xs- The standard deviation of 
these values was taken to be the measurement uncertainty 
in the mean observed On X s- We re f er to this value as the 
'bootstrap uncertainty' and denote it as Aboot(o"Nxs)- 

Estimates of the uncert ainty owing to the no ise pro- 
cess have been presented bv lVauehan et all i2003al) . based 
on light curve simulations, who showed that the noise uncer- 
tainty is proportional to the mean value of the variance. The 
constant of proportionality depends on the power spectrum 
shape, which we do no t know a priori. Therefore, as pointed 
out bv lVaughan et al 1 j2003aD . it is preferable to determine 
the uncertainties in o"nxs directly from the data. 

Our large database contains 6 objects (viz, AKN 564, 



MCG-6-30-15, TON S180, NGC 4151, NGC 3516, and 
NGC 5548) with a sufficient number of light curve seg- 
ments (>15) to make a meaningful estimate of the frac- 
tional uncertainty in cr^xs owing to noise nature of our 
light curves. First, we determined the standard deviation 
cr bs of the observed ct^xs values for each object. We then 
determined, from the bootstrap uncertainty, the standard 
deviation o" mea s that we would expect to observe in the dis- 
tribution of the crjsfxs values if the scatter was owing only to 
measurement uncertainties. We subtracted ff mcas from <7 bs, 
in quadrature, to obtain the standard deviation (T no i sc in the 
observed o"nxs values that is owing only to the stochastic na- 
ture of the variability. We then determined, for each of the 6 
distributions, the ratio between owse and mean cr^xs- We 
shall refer to this ratio as the 'fractional standard deviation' 
and denote it as af rac . The fractional standard deviations for 
the 6 objects were: 0.49 (AKN 564), 0.47 (MCG-6-30-15), 
0.69 (TON S180), 0.79 (NGC 4151), 0.82 (NGC 3516), and 
0.61 (NGC 5548). These values of erf rac show that, even in 
the absence of measurement uncertainties, one can expect 
noise uncertainties in the range ~50 -80 per cent for individ - 
ual measurments of o^xs ( see also lVauehan et al]l2003al) . 
This highlights the need of obtaining many realisations (i.e., 
many measurements of o"nxs)i regardless of the level of Pois- 
son noise in the data. 

IVauehan et"afl J2003al) showed that the uncertainty in 
the estimated variance of a red-noise light curve increases 
with the steepness of the power spectrum slope. Power- 
spectral analyses of AGN have revealed that the value of 
i'hfb generally decreases with increasing black hole mass. 
This means that the shape of the power spectrum in the 
frequency range probed by our light curves (~2.5 x 10 -5 
to 4 x 10 -3 Hz) is expected to vary as a function of M%, 
such that the objects with the highest M, should exhibit 
the steepest (a ~ 2) power spectra. We expect, then, 
that the scatter in cr^xs owing to red-noise fluctuations 
should also increase with mass. The lowest-mass objects are 
AKN 564 and MCG-6-30-15, and the observed values of 
ishfb for these fall within the frequency ran ge of our data 
iPapadakis et al.l2002l : IVaughan et al.l2003bl) . We therefore 
expect erfrac for this pair of objects to be less than the oth- 
ers. This does indeed appear to be the case: the crf rac values 
of AKN 564 and MCG-6-30-15 are both less than those of 
TON S180, NGC 4151, NGC 3516, and NGC 5548. However, 
we possess only a limited number of individual o^xs m ea- 
surements to estimate both the mean and standard deviation 
of each distribution, so it is possible that this apparent dif- 
ference is not statistically significant. We decided, therefore, 
to compare the 6 distributions of cr^xs values using a series 
of Kolmogorov-Smirnov (K-S) tests. Before we could com- 
pare the distributions, we first had to correct each of them 
for the effect of measurement uncertainties. To do this, we 
scaled the deviations of the observed cr^xs values so that 
the standard deviation of the 'corrected' distribution was 
equal to Cnoise, with the mean cr^xs remaining unchanged. 
We then normalised each of the corrected distributions by 
dividing the cr^xs values by the mean. The corrected and 
normalised cr^xs distributions for AKN 564 and MCG— 6- 
30-15 are consistent with being drawn from the same distri- 
bution. The same is true when comparing the other 4 dis- 
tributions with each other. We then created two, combined 
distributions: one for AKN 564 and MCG-6-30-15; and an- 
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other for TON S180, NGC 4151, NGC 3516, and NGC 5548. 
The fractional standard deviations from these two combined 
distributions were 0.48 (AKN 564, MCG-6-30-15) and 0.74 
(TON S180, NGC 4151, NGC 3516, NGC 5548), and a K-S 
test showed them to be different at the 95 per cent con- 
fidence level. The cumulative distribution functions of the 
combined distributions are presented in Fig. Combining 
the normalised, corrected distributions of all 6 objects re- 
sulted in a (Tf rac of 0.61. 

For the objects AKN 564, MCG-6-30-15, TON S180, 
NGC 4151, NGC 3516, and NGC 5548, we determined the 
total uncertainty [Atot ("nxs)] m the mean excess variance 
directly from their respective values of <7 bs- For each of 
the other objects, we estimated the noise uncertainty and 
combined it in quadrature with the bootstrap uncertainty 
[Aboot(oNxs)] using the following expression: 



Atot(o"NXs) — 



C fr. 



\ v 



+ [Aboot(of 



(2) 



where cr^xs i s the mean excess variance and N SC g is the 
number of available light curve segments. For objects with 
log M, > 6.54 we adopted a fractional standard deviation 
of ffrac — 0.74, while for the objects with log M% 5^ 6.54 we 
adopted a value of <7f rac = 0.48. These ranges in mass were 
selected on the basis that the object MRK 766, which has 
log M. = 6.54, is the most massive object that has an ob- 
served i^hfb in t he frequency ranged probed by our ~40 ks 
light curves (e .g., |P_apa£kkis_e±_alJl2Q^ fc Fabianl 

l2003l:lMarshall et alJl2004t IVaughan et alJl2004L and see~ 
troduction). For objects more massive than this we expect 
!/hfb to be less than our observed frequency range. In the 
absence of a measurement of M, or any information regard- 
ing the shape of the power spectrum, the mean value of 
ffrac = 0.61 can be adopted. 

The o"nxs upper limits for the non- variable objects were 
also estimated by combining the two components of uncer- 
tainty. We multiplied the la bootstrap uncertainty by the 
appropriate fractional standard deviation of the noise uncer- 
tainty. This value was then multiplied by 3 to provide an esti- 
mate of the 3cr upper-limit. We also estimated the 3cr 'boot- 
strap upper-limit' by multiplying the bootstap-uncertainty 
Aboot(o-Nxs) by 3. 

The distributions of the o"nxs values for AKN 564 and 
MCG— 6-30-15 are quite asymmetric, with each having an 
extended tail towards high values of (J^xs- However, the 
distributions of log erjvjxs look much more symmetric. This 
is not surprising as it is well known that the logarithmic 
transformation of a random variable with an extended tail in 
its distributi on brings that distribution m uch closer to 'nor- 
mality' (e.g.. lPapadakis fc Lawrencell993T) . Ideally, then, we 
would like to estimate log o^xs from each segment and then 
determine the mean of log onxs ^ or eacri object. Unfortu- 
nately, we cannot use this method because o"nxs i s negative 
for some light curve segments. We did, however, determine 
the logarithm of the mean onxs> which brings the distribu- 
tion of the mean cr^xs closer to normality. We also trans- 
formed the uncertainties A to t(o"Nxs) to be the uncertainty 
in the logarithm of the mean Onxs- 

The mean cr^xs values, uncertainties and upper limits 
are listed in Table The column listing Onxs gives the 



AKN 564 
MCG-6-30-15 




Normalised a.. 



Figure 1. 



Cumulative distribution functions of the combined 
distributions of AKN 564 and MCG-6-30-15 
(solid line) and TON S180, NGC 4151, NGC 3516, and NGC 5548 
(dahed-line). 



normalised cj^g 



uncertainty and 3cr upper limit as determined from only the 
bootstrap simulations. The uncertainties and upper-limits 
given in the column with log cr^xs include also the noise 
uncertainty. 



3.3 The variance— mass relation 

The relationship between log %xs and log M. is presented 
in Fig. |2| It is clear that there is a strong anti-correlation 
between the two quantities. This is confirmed using both a 
Spearman rank-order correlation test and Kendall's r, both 
of which show the anti-correlation to be significant with 
>99.99 per cent confidence. The upper limits to the vari- 
ance in the case where no variability is detected, which are 
shown in the upper panel of Fig. |5J are generally above the 
measured values (for a given mass) . This means they are un- 
likely to significantly affect any model fitting and we ignore 
them in the analysis below. 

While there is a strong general trend for objects with 
higher mass to be less variable, there is clearly substantial 
scatter in the variance-mass relationship. As we have, for the 
first time, presented realistic estimates of the uncertainties 
on o"nxs we can be confident that this scatter is not owing 
only to these uncertainties. 

There is also evidence from the plot-albeit based solely 
on the lowest mass object, NGC 4395-that the variance- 
mass relationship is non-linear. This is ex pected in the pres - 
ence of breaks in the power spectrum (e.g. . |Papadakisl2004h . 
as we now show by modelling the variance-mass relation- 
ship using both simple parametrizations and with a specific 
power-spectral form. 



4 MODELLING THE RELATIONSHIP 

BETWEEN EXCESS VARIANCE AND MASS 

Having obtained the mean cr^xs f° r each object, we then 
wished to model the relation between er^xs and M, . All fits 
were performed on log M, and log cr^xs- We fitted the data 
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using both a simple parametrization and also with a model 
that assumes the existence of a universal power spectrum. 



4.1 Simple parametrizations 

We initially modelled the data using a power-law of the form 
°nxs = AM^ 1 . The index and normalisation of the best- 
fiting power-law were 7 = 0.570 and A = 125, respectively, 
and the reduced chi-squared was xS/DOF = 8.05/31. This 
model is shown as a dot-dashed line in the top panel of Fig.|U 
We do not quote uncertainties in the best-fitting parameter 
values because the xi is formally unsatisfactory. 

We then used a singly-broken bending power-law de- 
fined as: 



°nxs = AM, 



1 + 



f— 

Km. 



m. 



Thigh — I'lo 



(3) 



where A is the normalisation factor and the function bends 
from a power-law slope of 7i ow to 7hi g h at the bend mass 

Af.bend- 

We fixed the lower index to 71 ow = 0. The best-fitting 
bend mass, normalisation, and upper index were M.bend = 
5.59 x 10 5 M Q , A — 0.144, and 7 high = 0.836, respectively 
(xl/~DOF = 5.99/30). The bending power-law clearly im- 
proves the fit statistic substantially, but it is difficult to as- 
sess the formal improvement with, e.g., an F-test as the fits 
are so poor. 



a functional form and is thus continous. We denote this con- 
tinuous power spectrum as Pm(v)- We need to determine 
how the discrete power spectrum Pd{v) is related to Pm(v)- 
The following description is appropriate for evenly sampled 
light curves containing no gaps and having an even number 
of bins. Note also that the model power spectrum Pm(v) 
must be defined to be two-sided and, since we are dealing 
with a noise process, we refer to the expectation value of 
each power. 

The first effect to consider is binning. Suppose we have 
a continuous process, with power spectrum Pm(y), and we 
transform it into a discrete process by binning the signal 
over a time period of St. The power spectrum of the ob- 
served binned light curve, say Pb(i'), is related to the Pyi(v) 
through the following relation: 

< Pa(y) >= B(y)Pu(y) (4) 

where the binning function B(v) Ivan der Klisll989l) is given 
by: 



B(y) = 



sm(jiv5i) 
TvuSt 



(•>) 



The next effect to consider is aliasing. The fact that the 
observed light curve is sampled at discrete intervals means 
that power can leak into the power spectrum from above the 
Nyquist frequency ^N yq = l/(2St). The binned and aliased 
power spectrum, say Pba.{v), is related to the intrinsic power 
spectrum Pm(v) through the relation «Priestlevl ll989 1 : 



2 

4.2 Predicting a NXS from a power spectrum model < p BA ^ >= < p B ^ + j/^) > 



(6) 



Based on recent power-spectral analyses of AGN, it is pos- 
sible that the power spectra of AGN have the same shape 
with the time-scale of the variations being proportional to 
black hole mass (see Introduction and references therein). 
We decided to investigate this possibility by modelling the 
relationship between cr^xs and M, with the assumption of 
a universal power spectrum. A model estimate can be made 
simply by integrating the continuous power spectrum over 
some frequency range, for example as defined by the length 
and the time bin size of the observation (e.g., IPapadakisl 
2004). This, however, neglects the effects from the sampling 
pattern of the light curve, specifically the fact that it is 
binned, may have gaps, and is of finite duration. For rea- 
sons discussed below these effects, particularly that of the 
finite duration and subsequent 'red-noise leak', are likely to 
be more important on the time-scal es considere d here than 
the much longer ones discussed by IPapadakisl i2004f) . We 
have taken an analytical approach to determining the model- 
predicted onxsi rather than us e simulations as is typical for 
power spectrum analysis (e.g.. ITlttlev et al"ll2002l) . Our ap- 
proach is preferable for two reasons. First, simulations are 
far more computer-intensive, and second they rely on the 
simulation technique accurately reproducing the character- 
istics of the physical process giving rise to the variability. 
While the technique described below applies to calculation 
of model cr^xs values it can be adapted straightforwardly to 
the estimation of discrete model power-spectra. 

According to Parseval's Theorem, the variance in a 
binned light curve is equal to the sum of the powers in 
the observed discrete power spectrum of that light curve. 
The model power spectrum, however, is initially defined in 



The power in one of our typical model power spectra de- 
creases sharply with frequency and the data are binned. 
This means that only a relatively small amount of power 
is aliased into the observed frequency range. Accordingly, 
we found that summing from i = — 10 to i = 10 was easily 
sufficient to account for aliasing. Power spectra that are ei- 
ther flat or increase with frequency might require a larger 
range in i. 

The final effect to account for is red-noise leak. This 
occurs when variations exist at frequencies lower than those 
sampled by the observed light curve, as is the case for a 
red-noise process. This 'leakage' of power from low to high 
frequencies can be seen as either a rising or falling trend over 
the duration of the light curve. The power spectrum of the 
final light curve, i.e. Pd(^), is related to the intrinsic power 
spectrum, i.e. Pm(v), by the convolution of the < Pba(v) > 
with the so-called 'window function' W(i>) of the observed 
light curve. For eve nly-sampled lig ht curves, W(y) is simply 
Fejer's kernel (e.g.. lPriestlevl ll989'l: 



W(y) 



1 
T 



sin(7ri/r) 



(7) 



where T is the duration of the light curve. We performed 
the convolution with the numerical integral: 

JV//2 

Pd(v) = 2 ^2 <Pba(zSv')> W{v-i5v')5v' (8) 

i=-Nf/2 

(v=l/T, 2/T,...,u Nyq ) 
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In the above sum, N is the number of bins in the light curve 
and / is a positive integer. The frequency step 8v' is given 
by 5v' = l/(Tf). The value of / must be large enough so 
that the convolution extends to a low-enough frequency to 
account for all of the low-frequency power. Determining a 
suitable value of / required a process of trial-and-error. We 
performed the convolution with successively higher values 
of / until further increases produced only a negligible effect. 
We found that / = 500 was sumcent for all our convolutions. 
Note that the introduction of the factor 2 in Eqn. |H| means 
that Paiy) is single-sided and it is defined only for N/2 
frequencies. Also note that, for i&u' = ±^N yq the term Sv' 
was replaced by Sv' /2, to account for the end-effects in the 
numerical integral. 

The expected excess variance was then determined by 
summing the powers in Pu(v): 



&NXS,model — 



'N/2-1 

Pu(i/T)Su 



(9) 



where 8v — 1/T. The factor of 1/2 is required for PD^Nyq) 
because in a double-sided power spectrum the Nyquist fre- 
quency occurs only once. The factor 8v is required because 
the power is expressed in units of fractional-rms-squared per 
Hz. 

Each of our 305 light curve segments has its own par- 
ticular duration and sampling pattern, and there are many 
gaps in the data train. Therefore, the window function will 
be different for each segment and will not, in general, be rep- 
resented by Fejer's kernel. However, the presence of missing 
bins in the light curve will affect only the scatter in the 
°nxs measurements, with the mean value being unaffected. 
Moreover, we have taken care to use light curve segments 
of similar durations. Therefore, we were able to simplify the 
modelling procedure by assuming that our light curves were 
all fully sampled with the same number of bins. We used 
N — 148, as this is the even number-of-bins closest to the 
mean segment duration of 38143.5 s. Having made this sim- 
plification, we were required to determine only a single value 
of 

°~nxs model f° r each object (for a certain model power 
spectrum), thus speeding up the modelling process. 

4.3 A universal power spectrum model 

Motivated by power- spectral analyses of A GN (see Intro- 
duction, in Darticulax jMarkgwit£_et al1 l2003). and following 
the recent work of |PapadakisTil2004[| . we hypothesised a uni- 
versal power spectrum of the form: 



Pyiiy) = A (^lfb/Vhfb) -1 ^ < ^lfb) 
Pm(v) = A (^/^hfb) _1 (^lfb < v < ^hfb) 
Pm(v) = A (^/^hfb)~ 2 (^hfb ^ v ) 



(10) 

(11) 
(12) 



where the normalisation factor A is the power at the high- 
frequency break ^-hfb- The value of i/hfb is assumed to 
decrease with black hole mass, according to the expres- 
sion 14jfb = Chfb/M,, where Chfb is a constant and M, 
is the mass of the black hole in units of Mq. The low- 
frequency break is related to the high-frequency break by 
z^lfb = ^fb/Clfb where Clfb is a constant. The normali- 
sation A varies as a function of ^hfb as A — PS Damp /ishfb , 
where PSDamp is assumed to be the same for all objects. 



Using this model, the relation between variance and mass 
can therefore be described with three parameters: Chfb, 
Clfb, and PSDamp- 

To determine the best-fitting model, we minimised x 2 
for grid of values of Chfb, Clfb, and PSDamp values. We 
found that we could not constrain the parameter Clfb- This 
is because the low-frequency break generally does not fall 
within our sampled frequency range. Therefore, we fixed 
this at Clfb = 20. This is rou ghly the value of Clfb ob- 
served in the AGN NGC 3783 iMarkowitz et al.ll2003D and 
in Cyg X-l in the l ow/hard state jBelloni fc HasingeJl990al: 
iNowak et aHll999h . 

The best-fitting values of Chfb and PSDamp are given 
in Table [5] This best-fitting model (for the fit including all 
33 objects) is shown as the solid line in Fig. (bottom). The 
probability of exceeding the Xv 01 t ne best-fitting universal 
model is 2 x 10 -25 . This indicates that, while the model 
appears to describe rather well the overall trend of decreas- 
ing cr^xg, there exists significant scatter not accounted for 
by the model. The residuals Alog onxs from this model 
are listed in Table We also fitted the universal model 
to the data excluding various objects. As seen in Table |5J 
neither the lowest mass object (viz, NGC 4395), nor the 
6 objects with the largest number of light curve segments 
(viz, AKN 564, MCG-6-30-15, TON S180, NGC 4151, 
NGC 3516, NGC 5548), dominate the fit. 

The scatter present in the relationship between log ctnxs 
and log M, can be explained with a variation of either Chfb 
or PSDamp from their best-fitting values. This is illustrated 
in Fig. [5] (bottom). We find that a range in Chfb values 
between 7.2 and 520 (upper and lower dotted-lines, respec- 
tively), or a range in PSDamp between 0.004 and 0.29 (up- 
per and lower dashed lines, respectively), can account for 
most of the scatter in the log cjsrxs versus log M, relation. 

The scatter might also be owing to a combination of 
the uncertainties in log cr^xs an d log M«, the latter of 
which are typically a bout 0.5 dex (e.g.. lWoo fc Urrvll2002l : 
iPeterson et alJl2004 - We performed simulations to inves- 
tigate this possiblity, adopting the best-fitting relation be- 
tween log o"nxs an d 1°S M m as our model. We needed first 
to obtain a set of 33 model data points to which we could 
then apply scatter in log cr^xs an d l°g M,. To do this, 
we projected each of our 33 observed data points onto the 
best-fitting relation, minimising the distance between the 
observed point and the model. (The distance between an ob- 
served data point and any particular location on the model 
relation was calculated from the differences in log cr^xs an d 
log Mb between the observed point and the model, divided 
by the corresponding uncertainty in the observed values.) 
Having thus adopted a set of 33 model data points, we then 
performed 1000 simulations. Each of these involved adding 
scatter to the model points and then determining the xi 
between the simulated data points and the model relation. 
We found that 79 per cent of the simulations produced a 
Xt exceeding that found for the observed data. Therefore, 
the scatter that we have observed in the relation between 
log cr^xs ana ~ m ighl De owing only to measurement 

uncertainties. If this is indeed the case, then we would ex- 
pect this scatter to be unrelated to other properties of the 
objects in our sample, and we investigate this possibility in 
the following Section. 
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Table 2. Best-fitting values for fits using the universal power spec- 
trum model. 



Excluded objects 
(1) 



Chfb PSD amp xJ/DOF 
(Hz M ) 

(2) (3) (4) 



None (all objects 43 
are included) 

NGC 4395 53 

AKN 564, 55 

MCG-6-30-15, 

TON S180, NGC 4151, 

NGC 3516, NGC 5548 



0.024 



i.24/31 



0.021 6.30/30 
0.033 4.30/24 



(1) Objects excluded from fit. (2) Scaling constant for the high- 
frequency break ^HFBi where j^hfb = Chfb/^«- (3) Power-spectral 
amplitude at i^hfb m power X frequency space. (4) Reduced chi- 
squared and degrees-of- freedom for fit. 




o tn 



4 6 8 

og M, (M Q ) 

Figure 2. Log of excess variance versus log of black hole mass. 
In the top panel, the dot-dashed and solid lines show the best- 
fitting power law and bending power-law models, respectively. In 
the bottom panel, the solid line shows best-fitting universal power 
spectrum model. The dotted and dashed lines illustrate the effect 
of varying either Chfb or PSDamp, respectively (see text for 
details). The o"^ xg upper limits are, for clarity, shown only in the 
upper panel. 
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Figure 3. Log of excess variance (top), log of the product of 
excess variance and black hole mass (middle), and excess variance 
residuals (bottom), versus log of the 2-10 keV luminosity. 



5 THE ORIGIN OF THE SCATTER IN THE 
VARIANCE-MASS RELATIONSHIP 

Previous studies have revealed an anti-correlation between 
°nxs an d X-ray luminosity, and a po sitive correlation be- 
tween Onxr and photon index V fe.g.. [Na ndra et alj|l997a|: 



[ Turner et a.Il 119991: iMarkowitz fc Edelsonll200lt iPapadakil 

2004). Given the strong dependence between the Onxs an d 
M, , it is of interest to see whether these correlations still ex- 
ist when this primary dependence is removed. This should 
allow us to shed light on the origin of the scatter in the 
variance-mass relationship. 

In Fig.[!J]we plot log o"nxs> l°g ^'^nxsi an< A the residu- 
als Alog o"nxs from the best-fitting universal model, versus 
the logarithm of the 2-10 keV luminosity. The quantities 
log M.a-^xs an d Alog On XS are useful because they remove 
the mass-dependence. Note that the quantity log M.ct^xs 
is model-independent. In Figs. 0] and [S] we plot the variabil- 
ity parameters versus, respectively, the photon index and 
the logarithm of the 2-10 keV luminosity normalised to the 
black hole mass, log (L2-10 ^ C y/M,). To the extent that the 
X-ray luminosity is proportional to the bolometric luminos- 
ity, as is commonly assumed, the value log (L2-10 ke\r/M m ) 
is proportional to the ratio between the mass-accretion rate 
and that required to reach the Eddington luminosity (i.e., 
the 'Eddington ratio'). Note that the correction factor be- 
tween the 2-10 keV and bolometric luminosities is uncertain, 
with considerable scatter. The Spearman rank-order corre- 
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Figure 4. Log of excess variance (top), log of the product of 
excess variance and black hole mass (middle), and excess variance 
residuals (bottom), versus the 2-10 keV photon index. 



lation coefficient and Kendall's r of all 9 relationships are 
presented in Table [2] 

As with previous studies, we find a very strong correla- 
tion between log cr^xs an d log Z/2- 10 k c v (see Fig.EJ. This 
correlation disappears when we remove the dependence of 
°nxs on M*. It seems most likely that the primary correla- 
tion is in fact with mass, and that the apparent correlation 
with log I/2-io kcv is secondary. 

A similar situation is present when considering the 
photon index (see Fig. 2J. Indeed, the correlation between 
log Onxs an d r is not very strong in any event, being signif- 
icant at only the 96 per cent confidence level, though there 
does seem to be an absence of objects having both a steep 
photon index and low (J^xs- When the mass dependence is 
accounted for, however, no residual correlation remains. In 
the plot of Alog o"nxs versus V, the steep spectrum objects 
do not have a systematically higher Alog o"nxs than the 
others. 

Finally, we consider the relationship between the 
variability properties and the normalised luminosity 
log (I/2-io kev/M.) (see Fig. 0. There is considerable 
scatter, and no strong correlation, between log cr^xs an ^ 
log (Z/2-10 kcv/M.). Here, however, we do find a signifi- 
cant relationship between log (L2-10 kcv/M,) and both the 
mass-normalised excess variance and the residuals from our 
best-fitting model. The latter correlation is significant with 
~99 per cent confidence and, perhaps surprisingly, it is in 
the sense that objects with larger values of normalised lumi- 
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Figure 5. Log of excess variance (top), log of the product of 
excess variance and black hole mass (middle), and excess vari- 
ance residuals (bottom), versus log of the 2-10 keV luminosity 
normalised by the black hole mass. 



nosity are less variable for a given mass. While significant, 
this relationship should be treated with some caution. The 
presence of random scatter in the black hole mass estimates 
could possibly induce such an anti-correlation. If M, is un- 
derestimated then Alog cr^xs wm a l so De underestimated 
and log (Z/2-10 kev/M,) will be overestimated. An artifi- 
cal anti-correlation would certainly be induced if all objects 
had the same value of log (Z/2-10 kcv/M.). However, it is 
less clear that this effect could produce an anti-correlation 
between Alog cr^xs an d (-^2-10 kev/M,) in our data 
because the normalised luminosities in our sample span 3 
orders-of-magnitude. We used the simulations described in 
Section 4.3 to test whether the observed anti-correlation 
could be owing to the uncertainties in the black hole masses. 
For each of the 1000 simulations, we calculated log cr^xs an< i 
log (Z/2_io kcv /M m ) from the simulated data points and mea- 
sured Kendall's r. We found that, even with no intrinsic anti- 
correlation between cr^xs an d £2-10 kev/M,, 57 per cent 
of the simulations gave a Kendall's r that was more nega- 
tive than the observed value of —0.31. Therefore, we cannot 
rule-out the possibility that the observed anti-correlation be- 
tween Alog erfSfxs an d l°g {L2-10 kev/M.) is an artifact in- 
duced by the presence of uncertainties in the measurements 
of black hole mass. 
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Table 3. Correlation coefficients between X-ray variability properties and the 2-10 keV luminos- 
ity, photon index and normalised luminosity. 





Observables 




Spearman 




Kendall 






Coeff. 


Sig. (per cent) 


Coeff. 


Sig. (per cent) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


log 


^2-10 kcV log cr^ xs 


-0.61 


99.98 


-0.43 


99.96 




log M.a NXS 


0.13 


53 


0.10 


59 




Alog 4 xs 


-0.06 


25 


-0.04 


26 


r 


log 


0.36 


96 


0.25 


96 




log M.cr^g 


0.10 


43 


0.10 


56 




Alog 4 xs 


0.11 


41 


0.09 


53 


log 


(^2-10 kev/M.) log er£ xs 


0.29 


89 


0.19 


89 




log M.4 XS 


-0.50 


99.7 


-0.36 


99.7 




Alog ^ xg 


-0.44 


99.0 


-0.31 


98.8 



(1) X-ray spectral property on the abscissa. (2) X-ray variability property on the ordinate. (3) 
Spearman rank-order correlation coefficient. (4) Significance of correlation. (5) Kendall's r. (6) 
Significance of correlation. 



6 DISCUSSION 

6.1 Summary of results 

We have investigated the relationship between normalised 
excess variance and black hole mass for a sample of 46 
radio-quiet AGNs. We restricted our light curves to have 
durations between ~30 and 40 ks (rest frame), allowing us 
to probe nearly the same range of time-scales for all ob- 
jects. There were 32 objects in our sample that had more 
than 1 light curve segment. For these objects, we were able 
to determine the mean cr^xsi decreasing the uncertainty in 
the measurements. Moreover, for 6 objects, there were more 
than 15 light curve segments available. An examination of 
the distributions of the individual onxs values for these 6 
objects allowed us to estimate the uncertainties in the mean 
°"nxs f° r every object in our sample. These uncertainties 
incorporate the effects of both measurement uncertainties 
and the stochastic nature of the variability. Of the 46 ob- 
jects in our sample, 33 were f ound to be var i able. As with 



previous studies using ASCA (Lu & Yu 2001; Bian & Zhao 


2003; 


iMarkowitz & EdelsorJ 12004!) and RXTE JPaDadakis 


El 


Markowitz & Ea^lson||2004|) data, we found a sienifi- 



cant anti-correlation between o"nxs and M m . 

We initially fitted the relationship between o^xs an< ^ 
M. with both a power-law and bending power-law. Neither 
of these fits were formally satisfactory, however the bending 
power-law was an improvement over the unbroken power- 
law. 

We also fitted the data with a universal power spec- 
trum model. We determined the expected o^xs from the 
model as a function of M. , accounting for the effects 
of binning, aliasing, and red-noise leak in the observed 
light curves. The best-fitting high-frequency-break x mass 
scaling-coefficent was Chfb = 43 Hz Mq, and the best- 
fitting amplitude wa s PSDamp = .024. In his study 
using RXTE data, iPapadakTsl i2004l) found values of 
Chfb = 1 7 and PSDamp = - 017 ( Chfb = 340 for 
NGC 4051). IMarkowitz fc EdelsorJ (120041) studied the vari- 
ability of Seyfert 1 galaxies on various time-scales and found 
that, on average, the variability time-scale followed the re- 



lation 7b = M./10 6 ' 7 days. Using our parametrization, this 
corresponds to a scaling factor of Chfb = 58 Hz Mq. 

In general, the mass- variance anti-correlation can there- 
fore be understood very simply by assuming that all size- 
scales scale with mass, and hence so do all characteristic 
time-scales (such as those represented by the break fre- 
quencies). Our analysis furthermore supports the idea that 
the average, or typical power spectrum of AGN resembles 
the 'universal' power spectrum discussed above. The best- 
fitting universal model was not satisfactory, however, with 
X^/DOF = 6.30/31, indicating that, for a certain M., there 
exists significant scatter in the o^xs values. However, our 
simulations showed that uncertainties in the mass measure- 
ments can account for this scatter. 



6.2 The origin of scatter in the variance— mass 
relation 

Previous work has suggested that the excess variance is 
related to source properties other than mass, such as 
the l uminosity, X-ray spectral index and H /3 line width 
(e.g.. lNandra et aljll997at iTurner et al.lll999l) . We have re- 
investigated some of these relations here. Consistent with 
previous work using ASCA data, we fo und a correlation be- 
tween log crjvxs and log L2-10 key (e.g.. lNandra et al.ll997al : 
ITurner et allil999l: I!7eighlvlll999ah . The fact that no corre- 
lation exists when the dependence of On X s on is re- 
moved suggests that the correlation between log Onxs an( A 
log I/2-10 kcV is largely a result of the o-^ xs -M, relation. 
This effect has also been seen in RX TE data with a time- 
scale of about 300 d (Paoadakis 2004). 

We also found an absence of objects having both a steep 
photon index and low o"nxs- After accounting for the de- 
pendence on mass, however, we found no evidence for a cor- 
relation between excess variance and X-ray spectral index. 
This is perhaps surprising, as previous work has suggested 
that narrow-line Seyfert 1 galaxies-whi ch have soft X-ray 
spectra as a gener al characteristic (e.g., iBoller et al] 119961 : 
iBrandt et allll997l) -are more variable than their broad-line 
analogues (ITurner et alll999tlLe7ghlvll999al) . An effect sim- 
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ilar to that which we have observed has already been noted 
by other workers using ASCA data. iLu fc Yul (1200 if) found 
that the narrow-line Seyfert 1 galaxies in their sample ap- 
peared to follow t he same varia n ce-ma ss relation as the 
broad-line objects. iBian fc Zhaol J2003I) . in an expanded 
study usin g the variance measurements of iTurner et alJ 
ll999l) and lLu fc Yul (I200ll1 . also found that the AGN with 
FWHM(RP) less than 2000 km s" 1 appeared to follow the 
same relation as those obj ects with broad H/3 emission lines 
(see also the discussion in iMarkowitz fc Edelsonll2004h . 

We also found an anti-correlation between ex- 
cess variance residuals and the normalised luminosity 
(L2-10 kcv/Af.), which we shall now simply refer to as the 
Eddington ratio M. Our simulations showed that this ap- 
parent anti-correlation between Alog o^xs an d M could be 
an artifact owing to the uncertainties in the measurements 
of the black hole masses. The fact that we did not find a 
positive correlation between excess variance and M, for a 
given mass, is surprising: in the prevailing paradigm, NLSls 
generally show more variability and are thought also to be 
accreting at high Eddington ratios iPounds et alJll995ft . It 
is not yet clear, then, that a high value of M is a contribut- 
ing factor to an AGN exhibiting a relatively large excess 
variance. Further investigations in this regard will benefit 
enormously from future improvements in black hole mass 
measurements. 

6.3 Models for X-ray variability 

In the standard coronal model, which can be applied 
both to stellar-mass black holes and AGNs, seed photons 
from an optically thick accretion disc are inverse Comp- 
ton s cattered by hot electrons in an accretion disc corona 
fe.g..lSunvaev fc Titarchuklll980l lHaardt fc Maraschill993t 
IChurazov et alJl200lt iMcClintock fc Remillardll2004l) . 

One class of models involve s the superpo sition of indi- 
vidual 'shots' in the light curve llTerrel]lll972h . These shots 
arc possibl y associated with m agnetic flares in the corona 
fe.g.. IPoutanen fc _Jiabian|LL999l and refe rences therein). In 
the model of IPoutanen fc Fabianl IJ1999I) . there is a distri- 
bution of shot time-scales, with the value of ^hfb being in- 
versely proportional to the duration of the longest shots. 
Also in that model, the variance of the counting rate fluctu- 
ations is inversely proportional to the mean rate A of the oc- 
currence of flares. One can then assume a basic framework in 
which all size-scales (and, therefore, time-scales) and the lu- 
minosity of the individual shots is proportional to the black 
hole mass, accounting for the main variance-mass relation- 
ship. The total luminosity is proportional to A, so for a given 
black hole mass the variance in the light curve is expected 
to be inversely proportional to the Eddington ratio. 

In the so-called 'propagating pertubation' class of 
models, variations in the accretion ra te occur over a 
range of radii from the black hole (e.g. , |Lvjibarakii[ Jl997l: 
IChurazov et al J 1200 it iKotov et all l200lt lUttlevI 12004. and 
references therein). Slower variations occur at larger radii 
and propagate inwards, coupling together with the faster 
variations produced at smaller radii. The modulations in the 
accretion rate propagate to the X-ray emission region and 
produce variations in the X-ray flux. This type of model is 
attractive because it can provide an explanation for the well- 
know 'rms-flux' relation seen in X-ray binaries and AGN 



(e.g., lUttlev fc McHardvll200ll Illttlevll2004t iGaskeljEooil . 

The value of vhfb is expected to be inversely proportional 
to the size of the X-ray emission region because the varia- 
tions th at originate from within the emissio n region are sup- 
pressed Jchurazov et al.ll200ll:ITittlevl2004l) . In the model of 
IChurazo^Tt^iTl20(nF ~the low/hard state in Cyg X-l occurs 
when the optically thick, geometrically thin accretion disc 
is truncated far from the emission region. In the high/soft 
state, the disc reaches all the way down to the emission 
region and this leads to the X-ray variations following an 
unbroken a = 1 power-law. In this model, it is not fully 
specified how the emission region changes as the inner ra- 
dius of the disc varies. It is clear, however, that the emission 
region would need to become smaller as the disc approaches 
that region because vbfb is higher in the high/soft state 
tha n in the low/hard state . 

iMcHardv et alJ <l2004h appealed to the analogy with 
black hole X-ray binaries and speculated that the location 
of the inner edge of the accretion disc in AGN is perhaps 
related to the mass-accretion rate or the black hole spin. 
For a certain black hole mass, then, different AGN might be 
regarded as existing in different states, just as Cyg X-l is ob- 
served in different states. In this scenario, we would expect 
the X-ray variability of AGN to be related not only to the 
black hole mass but also the Eddington ratio and photon 
index. Objects having a relatively high M and soft X-ray 
spectra would, for a certain value of M,, have a relatively 
high value of vhfb (i.e., a high value of Chfb) and should, 
therefore, exhibit a relatively high value of o"nxs f° r a given 
range in time-scales. We found no evidence that the X-ray 
variability depends on these properties, and so the reality of 
this scenario remains to be established. Note, however, that 
if an ant i- correlation existed between Chfb and PSDamp, 
then Chfb could possibly increase without there being a 
corresponding increase in cr^xs- 

Discriminating between various possible scenarios obvi- 
ously requires the use of power spectral analyses, preferably 
covering a wide range in source properties. The challenge, 
then, is to assemble enough high-quality power spectra so 
that we can relate the power-spectral parameters not only 
to M. but also to Eddington ratio and other quantities such 
as photon index. We note, in particular, that an analysis of 
the AGN data in the XMM-Newton and Chandra archives, 
even from relatively short observations, would be useful in 
studying the properties (e.g., power-law slopes) of the vari- 
ability at frequencies above the high-frequency break. A nat- 
ural starting point, of course, is to conduct a rigorous com- 
parison between the currently available power spectra (e.g., 
l Uttlev et all 120021: IMarkowitz et all 120031: IMcHardv et alJ 
2004) and the other relevent source properties. Any conclu- 
sions draw from these comparisons could then be tested on 
a larger sample of objects by using measurements of excess 
variance. 
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